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Abstract 

We prove existence of global in time strong solutions to the truncated thermo- 
visco-plasticity with an inelastic constitutive function of Norton-Hoff type. This 
result is a starting point to obtain renoramlised solutions for the considered model 
without truncations. The method of our proof is based on Yosida approximation of 
the maximal monotone term and a passage to the limit. 


1 Formulation of the problem and main result 

Our study is directed to mathematical analysis of thermo-visco-plasticity (for derivation 
we refer to [22, 24] and ESI). This means to problems from the theory of inelastic de¬ 
formations in which the temperature affects the visco-plastic response of the considered 
material. Let us assume that the elastic constitutive stress-strain relation has the form 

a = C(e(u)-e p )~fm, (1.1) 

where o is the Cauchy stress tensor, e(u) = |(Vu + V T u ) is the linearized strain tensor, 
u is the displacement vector, e v is the inelastic strain tensor, 6 is the temperature, f is a 
given function depending on the considered material and C is the elasticity tensor which 
we assume to be symmetric and positive definite on the space of symmetric matrices. 
Notice that the thermal part of the stress —/(0)11 is not linearized in the neighborhood 
of the reference temperature (compare [5]-[8] and [291150] ). Our motivation for the form 
of the elastic constitutive relation (II. ip follows the results of [12j| and [13] . 

Moreover, we assume that the density of the internal energy e has the simple form 
e = c9 + C ~ l T ■ T, where T = a + f(9) U. In this case as a consequence of the first 
principle of thermodynamics we obtain the following form of the heat equation 

cp9 t — xA9 + f(9)div u t = T ■ e v t , (1.2) 
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where p > 0 is the mass density, k > 0 is the material’s conductivity. It is also assumed 
that the evolution of the inelastic strain tensor e p is given in the form 


4 = G(T), 


(1.3) 


where G is a given maximal monotone vector held with G(0) = 0 (inelastic constitutive 
equation with only one internal variable £ p , see In the main part of the article we 
will specify G choosing it in a form of the Norton-Hoff model (similar nonlinear how rule 
was considered in |3j). 

If we consider equation (11.21) with the homogeneous Neumann boundary condition and 
the homogeneous balance of forces div a = 0 also with homogeneous boundary conditions 
then we can conclude that the considering problem possesses a natural semi-invariant 
function namely the total energy 


£{t) = 


cp9dx -\— / C l T ■ T dx < £(0), 


(1.4) 


where h C M 3 is a bounded domain describing our considered body with boundary of 
class C 2 . In this simple calculation we have used the dissipation inequality, which yields 
that T ■ > 0 in the whole deformation process. From the observation (jl .4[) we see that 

the temperature is controlled in the space L 1 (h2) only (with additional information that 
6 > 0). Moreover, in general the term T ■ e v t belongs to L 1 (fl) only. From these reasons 
we are going to prove existence of solution in renormalised sense -see for example 0-m 
and [13, IH [16, ,211 EE]. In the literature there are only some articles with mathematical 
study of special thermo-visco-plastic models with various modifications (see for example 
[5]- [7] andpjH 23]). See also the articles [57J and [28], where a poroplasticity models are 
investigated which have a similar structure to the linear thermo-plasticity. 

In this article and in the following work [[2(0], we are going to study the therm-visco¬ 
plasticity model of the Norton-Hoff type with a damping term, which we interpret as 
external forces acting on the material and depending on the deformation velocity. Thus 
the system of equations, which we study in this article is in the form 


diva, a 
a 

F P 

T 


9 t - Ad + f(9)divu t 


-F - div C(e(u t )), 

C (s(u)-e p )-fm, 

| dev(T) | p_1 dev(T), 
C(e(u)-e*), 

| dev(T)| p+1 , 


(1.5) 


where p > 1 is a given real number and dev (T) = T — | tr (T) • 11 is a deviatoric part of 
the symmetric matrix T. The function F : H x [0, T] —» M 3 describes the density of the 
applied body forces. The main idea in the existence theory of renormalised solutions is 
to study the so called truncated problem and next to prove that the sequence of obtained 
solutions converges to a renormalised solution. 

In this article we study the truncated model only, hence for fixed T > 0 and e > 0 we 
have to find the displacement held u : fl x [0, T] —> M 3 , the temperature of the material 
9 : Q. x [0, T] —> M and the visco-plastic strain tensor e p : x [0, T] —> S^ ev (5(j cv denotes 
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the set of symmetric 3 x 3-matrices with vanishing trace) satisfying the following system 
of equations 


cliv x a = -F - cliv C(e(u t )), 

a = C(e(u)-£n-f(T±(m, 

s? = | dev(T) \ p ~ 3 dev(T), (1.6) 


T = C (£{u)-£ p ), 

6 t -Ad + f(Ti_(9))divu t = Ti(|dev(T)| p+1 ), 

where the function Ti (•) is the truncation at height - > 0 i.e. Ti (r) = min(i, max(r, —-)). 
In this paper we assume that the function / £ C' 1 (M; M). The system (11.6(1 is considered 
with the Dirichlet boundary condition for the displacement and the Neumann boundary 
condition for the temperature 


u(x,t ) = g£>(x,t) for x £ and 1>0, 

d 9 

-Q—(x,t) = go(x,t) for x £ Oil and 1 > 0. (1.7) 

Finally, we consider the system (1 1.6 [1 with the following initial conditions 

u(x, 0) = Uo(x), £ p (x, 0) = £ p,0 (x), 9(x, 0) — 9 0 (x). (1.8) 

Suppose that our data have the following regularity 

F e H 1 (0,T;L 2 (fl;R 3 )), g D £ bF 2 ’°°(0, T; (0fi; M 3 )), (1.9) 

geeH\0,T-L 2 (dn-,R)), (1.10) 

u 0 £ H\£l-R 3 ), £ p ’° £ L 2 (fl; S 3 ev ), do £ H 1 ^, R) . (1.11) 

Next, we define a notion of a solution for the system (11.61) . 


Definition 1.1 Suppose that the given data satisfy \1.9d - \l.ll\) . We say that a vector 
u £ if 1 (0,T; H^Q^R 3 )), the function 9 £ L°°( 0,T; H\Q;R)) such that9 t £ L 2 (0,T; L 2 (fl;M)) 
and the inelastic deformation tensor £ p £ H 1 { 0,T; L 2 (Q;<S| ev )) are solutions of the prob¬ 
lem U.60 - U.81 ) if the equations (ll.6jl i - (jl.6(1 /i are satisfied for almost all (x, 1) £ 12 x (0, T ) 
and the heat equation is satisfied in the following sense 

f 9 t v dx + f 'V9'Vvdx + / f(Ti(9))divu t vdx = f Ti (| dev(T)| p+1 ) vdx 

Jn Jn Jn e Jn e 

+ / g 9 vdx 

Jan 

for all v £ H 1 ( 12; R) and almost all t £ (0, T). 


Theorem 1.2 (Main Theorem) 

Suppose that the given data have the regularity required in l \1.9\) - in~Tl\) . Moreover, let 

| dev(C(e(M 0 ) - e p ’°)| p " 1 dev(C(e(u 0 ) - £ p ’°) £ L 2 (Q; S$ ev ). 

Then, for all T > 0 the system 41.61) with the boundary condition and the initial 

condition 111.81) possesses a solution in the sense of Definition 11.11 
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The initial assumption in Theorem 11.21 means that for t = 0 the argument of the maximal 
monotone operator belongs to the domain of | dev(T)| p_1 dev(T) (for more information 
we refer to [2], Definition 2.4 and Theorem 2.5). 

Our Theorem 11.21 will be proved in the next three sections. First, we use the Yosida 
approximation to the maximal monotone inelastic constitutive equation. Then we pass 
to the limit to obtain a solution in the sense of Definition 11.11 


2 Transformation to a homogeneous boundary-value 
problem with respect to the temperature 


Let us consider the following linear parabolic system 

9 t (x, t) — A9(x, t) — 0 (2.1) 

for x G dCl and t > 0 , 

for ieO. (2.2) 

Assuming that gg G iY 1 (0, T; L 2 (<9D; M)) and 9 o G iL 1 (0;M) we conclude that the sys¬ 
tem (12.11) possesses a solution 9 G L°°(0, T; iL 1 (D; M)) such that 9 t G L 2 (0, T; L 2 (D; R)). 
Additionally the following estimate 

II || L 2 (0,T \L 2 (fi;R)) + || 9 1| L°° (O.T;# 1 (0#)) < D ( || Qg || H 1 (0,T;L 2 (5O;R)) + || 9q || H 1 (Q;R) ) (2.3) 


with boundary-initial conditions 
89 


(x,t) = g 0 (x,t) 
8n 

9(x,0) = 9 0 (x) 


holds. Finally, we define 9 = 9 — 9. Notice that to find a solution (u, 9) to the problem 
(II.6p ~ (II. 81) we need to find a solution (it, 9) of the following problem 


div x a 
o 

F P 

t t 

T 

9 t - A9 + f(Ti(9 + <9))divu t 
with the initial-boundary conditions 

u \an 9 d j 


0(0) — 0o — i 


= —F — divC(e(it t )), 

= C(e(it) — e p ) — f (Ti (0 + 0)) 11, 

= | dev(T)| p_1 dev(T), 

= C (s(u)-s p ), 

= Ti(|dev(T)| p+1 ) 

e 

89 

IT =0 ’ 

8n |an 

lt(0) = U 0 , e P (0) = £ p ’° . 


(2.4) 


(2.5) 


3 Yosida approximation 

Notice that the inelastic constitutive equation (12.4ft ^ is a maximal monotone vector held 
such that G'(O) = 0, where G(T ) = | dev(T)| p_1 dev(T). Moreover, there exists a positive, 
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differentiable convex function M : 5| ev —> <S| ev such that S7 T M = G. The main idea 
to prove an existence result for system (12.4p is based on the, so called, partial Yosida 
approximation. This means that we will use the Yosida approximation of the maximal 
monotone vector held VM. For A > 0 let us define the function 

M\(z) = inf -^-{\z - w\ 2 + M(w)}. 
toeR 6 2A 


Hence, M\ is subquadratic, nonnegative, the gradient VM A is a global Lipschitz map and 
is the Yosida approximation of VM - for more information we refer to [4] and H3- Using 
the function M A we define a sequence of approximate problems 


div x a x = — F — div C (£ X (u t )), 

a x = C{e{u x )-e p ’ x ) - f(Ti_(6 x + 

£ P t X = VM A (dev(T A )), 

T x = C(e(n A ) - £ p ’ x ) , 

9 X - A9 X + f(Ti(9 x +~6)) div u x = Ti(dev(T A ) £ p,x ) . 

The problem (13.ip is considered with the boundary conditions 

u x (x,t ) = go{x,t) x G dfl, t > 0, 

de x 


dn 


(x,t) = 0 


x E cKl t > 0 


and with the initial conditions 

u x (x, 0) = uq(x), £ p ’ x (x, 0) = e p ’°(x), 0 X (x, 0) = d Q (x). 


(3.1) 


(3.2) 


(3.3) 


To prove existence result for all A > 0 we are going to define the operator 1Z acting from 
the Banach space L 2 (0, T; ih 1 (h2; M)) into the same space. Next, using to this operator 
the Banach Fixed Point Theorem we find a solution of system (13.ip with initial-boundary 
conditions (13.2p and (13.311 . 

Fix A > 0. Let us set 9* E L 2 (0, T; iL 1 (A2; M)) and consider the visco-elasticity problem: 

div^cr = — F — div C(e(tq)) , 

a = C(e(u) -£ p ) - f(Ti(9* + 6»))l, 

£ p = VM A (dev(T)), (3.4) 

T = C{£(u)-£ p ), 

u \en 9 d i 
u( 0) = u 0 , 
e p (0) = £ p ’°. 

Notice that in the system (13.4p we dropped the subscript A > 0. In the next part of this 
section, we will drop the subscript A and write u, 9 and £ p instead of u x , 9 X and £ P)X . 
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Lemma 3.1 Let us assume that the given data satisfy all requirements of Theorem \1.2\ 
and 6 * G L 2 (0,T; Then, there exists a global in time unique solution (u,£ p ) of 

the system ffFff such that 

(u,e p ) G HfO^Hf^R 3 )) x T; L 2 (F>; S 3 ev ))- 

Proof: The idea of the proof is a fixed-point argument. Let £* G L°°(0,T; L 2 (Ll]S 3 )) 

and let us consider the equation 

e p (x, t) = £ p,0 (x) + J VMx^dev (C(e*(r) — £ p (r))) jdr . (3.5) 

From the theory of differential equations in Banach spaces (VM\ is global Lipschitz) it 
follows that the equation (13.5ft possesses a global in time, unique solution 
£ p G W 1,oo (0,T; L 2 (Q- «Sj ev )). Next, consider the system of equations 

div^cr = -F - divC(e(tq)) , 

a = C(e(u) -£ p ) -/(Ti(0* + 0))lL, (3.6) 

u(0) = u 0 , u\ gn = g D , 

where £ p satisfies (13.5p . Notice that the term f(Ti(9* + 0)) is bounded, because the 
argument of the continuous function / is bounded. Therefore the system (j3.6[) has a 
unique solution u G H 1 { 0, T; H l {Vt] M 3 )). This way, we defined an operator 

P : L°°(0, T; L 2 (fi; S 3 )) L°°(0, T; L 2 (Ll; S 3 )) (3.7) 

such that P(e*) = |(V w + V T w). We will show that P is a contraction. 

Let us denote by £i(t) and £%(t) solutions of (13.5ft with the input functions 
£*, £* 2 G L°°(0,T; L 2 (Q]S 3 )), respectively. The difference P(£*) — P(£ 2 ) satisfies 

- div (C(P(£t) - P(£* 2 ) + P(£l) t - P(£*) t )) = -div (C(£? - 4)) (3.8) 

and the following inequality 

Ili’KW) - ^KWJIUw) < C|K(<) - 4MllLWi.) (3.9) 

holds, where the positive constant C does not depend on these input functions, 6 * and is 
independent of t. Using equation (13.5p it is not difficult to obtain the following inequality 


I4W - ejWIUwij £ dkiW - 4( ( )IUw»)’ 


(3.10) 


where C does not depend on t (it depends only on the Lipschitz constant and on time 
T). Having these two inequalities we easily see that the operator P is a contraction- for 
more details we refer to [2].. 

□ 

Theorem 3.2 Suppose that the given data satisfy all requirements of Theorem \1.2[ For 
all A > 0 the system \3.1\) with initial-boundary conditions \S.‘d\) and hS.tft) possesses 
unique, global in time solution ( u, £ p , 9) such that 

(u, £ p ) G H\0, T; M 3 )) x W^O, T; L 2 (H; 5 3 ev )), 

9 G L°°(0, T- M)) and 9 t G L 2 (0, T; L 2 (fi;! 
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Proof: Fix 9* G L 2 (0, T; H l {Vt] M)). Lemma 13.11 implies that there exists solution 

(w, £ p ) G LU(0, T; M 3 )) x hF 1,oo (0, T; L 2 (fi; <S| ev )) of the system (13.4p . Inserting this 

solution to heat equation (I3.1|L we obtain a solution 9 G L°°(0, T; M)) such that 

G L 2 (0,T; L 2 (0; M)). Hence, we have defined an operator 


L 2 (0, T; M)) 7£(0*) = 9 G L 2 (0, T; 

Next we prove that 1Z is a contraction. Let (iq, ef) and (« 2 , e 2 ) be the solutions to (13. 4 p for 
9\ and 9\ G L 2 (0, T; iL 1 (Q; R)), respectively. Let us mark by u = iq — u 2 and e v = £ p — £ p . 
Then £ p (0) = 0 and h(0) = 0. The weak formulation of the system (13.4p yields 

f C(£(u) + £(u t ) — £ p )e(v) dx = f (^f (T i (9 1 + 9)) — f(T i (9 % + 9))^J div v dx (3.11) 

for all v G Lfo(f2,M 3 ). Putting v = u t in (13.lip we obtain 
1 d 


2 dt 


in 


C(e(u) — £ p )(e(u) — £ p ) dxj + / C £{u t )£(u t ) dx = — C(e(-u) — £ p )£ p dx 


f{T±(6\ + 9)) - /(Ti(0* + Oj) )divu t dx 


(3.12) 


Notice that the difference f(Ti(9^ + 9 )) — f(T\ [9^ + 9)) is bounded. Using monotonicity 
of VMa and integrating with respect to time we get 


C £(u t )£(u t ) dxdt < T C 


(3.13) 


'o Jn 


and the constant C > 0 does not depend on 9\ and 9\ (it depends on e and H only). 
Additionally, from the weak formulation of (13 .1 p we have 


1 d 

2 dt 


'n 


( J \9, - 9 2 \ 2 dx) + J \V(9 1 - 9 2 ) | 2 dx = 

f{T±(9 x + 0))div u\ - f(Ti_{9 2 + 9))di vw 2 )^ - 0 2 ) dx 

Ti (dev(Ti) sf) - Ti (dev(T 2 ) £ p ’ 2 )) ( 9, - 9 2 )dx , 


(3.14) 


where 9 1 and 9 2 are solutions of (13.ip ^ with the input functions (ui,£ p ) and (u 2l £ p ), 
respectively. Using the Cauchy-Schwarz inequality with small weight and integrating 
with respect to time we get 


\9\ — 02 1 2 dx + 


|V(#i — 0 2 )| 2 dxdt < C(v) 


| £{u\ — u 2 )\ 2 dxdt 


in 


i o Jn 


I o Jn 


|6>! — 6> 2 | 2 dxd t + TC, 


(3.15) 


'0 Jn 


where the constant C > 0 does not depend on 9\ and # 2 . Inserting (I3.13P into (I3.15P and 
choosing v > 0 sufficient small we receive 

[ |6>i — 6> 2 | 2 dx + [ T [ |V {9\ — 9 2 )\ 2 dxdt <T D f \\9\-9l\\ 2 mm df (3.16) 
Jn Jo Jn Jo 
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where the constant D > 0 does not depend on 9\ and 9\ (it depends on e and only). 
Moreover, notice that D is also independent on the initial data e p, °, u$ and 9o, thus the 
inequality (j3. 16|) implies that the operator 1Z is a contraction. □ 


4 Passing to the limit in Yosida approximation 

This section is devoted to prove some a priori estimates for the sequence of approximate 
solutions {(u A , 9 X , £ p,a )}a>o and pass to the limit A —> 0 + . In this section we returned to 
the subscription A. 

Remark: The construction of the operator 1Z yields that 9f G L 2 (0,T; L 2 (A2; R)) and 

we can conclude that the following elliptic problem 

div(c( £ (« A ) - £p) - /'(Ti(0 A + «))vr t (» A + g)(6t + a,)l) = —F t - divC(£( M A )), 

u t In = dD,t 

possesses the solution u x G H 1 ( 0, T; R 1 (0; R 3 )) (Ti is a Lipschitz function and the chain 

e 

rule for the derivation of Ti holds true). 

e 

Let us start with prove some estimates for the approximation sequence. 

Theorem 4.1 (energy estimate) 

Assume that the given data satisfy all requirements of Theorem \1.2[ Then there exists a 
positive constant C(T ) (not depending on X) such that the following inequality holds 


\\T\t)\\ 2 L2{n . s3) + / ||e A (r)||i 2(n;<s3) dr < C(T) 

Jo 

for t <T. 

Proof: Calculate the time derivative 

j (i f C(e A - e p,A )(e A - e p - A ) dx) = / C(e A - e p ’ A )(e t A - £ p t ’ x ) dx 

r ' A e A dx+ [ f (Ti (9 x + 9)) div u A dx — [ dev(T A ) ef A dx . 


a 


(4.1) 


The last integral on the right-hand side of (14.ip is non-negative. Integrating with respect 
to time we obtain 


C(e A - e p,A )(e A - e p ’ A ) dx < / C(e A (0) - e p ’°)(e A (0) - e p ’°) dx 


'n 


in 


+ 


r a x e A dx dr + [ [ f(Ti {9 X + 0)) divn A dxdr. (4.2) 

>o Jn Jo Jn 
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Integrating by parts in the second term on the right hand-side of (j4.2[) and using the 
equation (13.lib we get 



o Jn 


/ C(e A -e p ’ A )(e A -£ p ’ A )dx + 

Jn 

< [ C(e A (0) - e p ’°)(£ A (0) - £ p ’°) dx + 

+ C(£(u A )))ng D)t dS dr + 


C(e A )e A dx dr 



F u+ dx dr 


0 J n 


+ 



o Jan 


a 



o J n 


f (Ti (9 x + 0)) div u x dx dr . (4.3) 


The boundary integral appearing on the right-hand side of inequality (14. 3 p can be esti¬ 
mated using the continuity of the trace operator in the space L 2 (div) [27] 


(a A + C {e{u x )))ng Dtt d S dr 



o Jon 


< c 


a 


+ <C(e(u t ^ r ) 


H U (3Q;R 3 ) 


< c( / (||a A + C(e(u A ))|| L 2 (n . S 3 ) + ||div(a A + C(e(u A ))|| L 2 (n . R 3))||c/ A 


4,1 .Hr (30;R 3 ) ^ r 


— C ( / (ll°" A + ^-'( e (' U t'))lU 2 (n ; 5 3 ) + H-^l|L 2 (0;R 3 ))ll5 , D,tlL,dr 


'H?(d£l; R 3 ) 


< C 


+ 


T 11x2(0-53) + ||i ? || i 2(o ; R3)) ||g£> )t 


k A |U 2 (0;5 3 )lbD,t|| H i (an;K3) dr + 1 


i dr 

R 3 ) 


(4.4) 


Putting (1 1. Hi into (14. 3 p we receive 


C(e A — e p,A )(e A — £ p,x ) dx + / / C(e A )e A dxdr < C(T) 


o Jn 


< / C(e(u 0 ) — £ p, 0 )(£(u 0 ) — £ p ’°) dx + v / / |e A | 2 dxdr + v \T X \ 2 dx dr , (4.5) 


o Jn 



o Jn 


where v > 0 is any positive constant and C(T) does not depend on A. Finally, choosing 
in (14. 5 p v > 0 sufficient small we finish the proof. □ 

Remark: The Theorem 14.11 implies that the sequence {m a }a>o is bounded in the space 

L 2 (0,T; i7 1 (f2;M 3 )). Let us consider the equation (13. ip ^ 

9 X -A9 X = Ti(dev(T A )ef A ) - f(Ti(9 x + 9))dwu x . (4.6) 

Notice that the right hand-side of (14.6p is bounded independently on A in the space 
L 2 (0, T; L 2 (f2; K)). From the standard theory for heat equation we immediately conclude 
that the sequences {# a }a>o an d {# A }a>o are bounded in the space L°°(0, T; H 1 ( fl; M)) and 
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L 2 (0, T; L 2 (fl; K)), respectively. 


The next step in onr existence theory is an estimate for time derivatives of the approximate 
sequence. The proof of the theorem below works for gradient flows only. The main idea 
is similar to that from Theorem 2.3 in [TBj. 

Theorem 4.2 (L 2 (L 2 ) estimates for time derivatives) 

Assume that the given data satisfy all requirements of Theorem \1 .A Then, for all T > 
0 and t < T the solution of approximate problem hS.l\) - \3.‘A) satisfies the following 
inequality 

t 

j H^ A ( r )llL 2 (t2;5 3 ) d r + ll £A ( r )llL2(Q ;< S3) < C(T), (4.7) 

0 


where the constant C(T) does not depend on X. 
Proof: Let us fix T > 0. Compute the derivative 
d 
dt 


M A (dev(T A ))dx = / VM A (dev(T A ))(dev(T A ))da; 


dev(C( £ K A )-ef A )) £ f A dx= / T A sf A dx = - / C^T^dx 


in 


+ 


T t £ t = 


C 1 T, X T. X d.r + / erf ef dx 


A ,-A. 


t t 


t 


+ / ./'(vi w' + »i)vr, (» A - ») y + »,jfiiv..;dr. 


(4.8) 


Notice that in the formula (14.81) we use the following informations: Ti is a Lipschitz 

e 

function, Ti{6 x + 6 ) belongs to 77 1 (h2) and the chain rule for derivation of Ti(0 x + 6 ) 
holds true (see e.g. m Esi). Integrating the equality (14.8ft with respect to time, we get 


M A (dev(T A (t)))dx + 


in 



o Jn 


C- 1 T A T A dxdr+- / C(ef)sfdx 


M A (dev(T A (0))) dx + - / C(e A (0))e A (0) dx + 

^ ./O 



o Jn 



(y + C(4))e?<MT 

(4.9) 


o Jn 


f(Ti(d x + 6))VTi_ ( 9 X + 6) (9f + 9 t ) div uf dx . 

The function rq(0) is the solution of the following linear elliptic problem 

divC(e(n A (0))) = -F(0)-div(c( £ M-£ p ’°)-/(Ti(d 0 ))ll 

u x (x, 0) = u 0 (x), «t(0)| 8n = g D ,t( 0) • 

Therefore, we obtain inequality 

ll U t (0) II 1/ 1 (£7;M 3 ) — ^(ll^ollHqOiR 3 ) + IIIIL 2 (0;R 3 ) + 11^’°IIL 2 (0;S| ev ) 

2 

H?(dn; R 3 ) 


(4.10) 


\\f{Tl(0o))\\l H n m + \\9D, t m 


(4.11) 
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The inequality (14.lip implies that ||u A (0)||iji(n ; R3) is bounded (independently of A). The 
function Ti is Lipschitz, hence VTi is bounded. Moreover, the remark before the Theorem 

e e 

EO implies that the last term on the right hand-side of (14. 9p is bounded. Now we integrate 
by parts in the third term on the right hand-side of (14.9p to get 



o Jn 


[a x + C(e A t )) e x dxdr = / / F t u x dxdr + 


o J n 



o Jan 


(a? + C (e^))ng D>t dS dr 



F t u} dx dr 


+ [ + < £(4( t ))) n gDA t ) dS - [ (cr A (0) + C(e t A (0)))n^ Di i(0) dS. (4.12) 

Jan Jan 

The boundary integrals are estimated using the continuity of the trace operator in the 
space L 2 (div) (in the same way as in (14.41) 1 hence we get 

[ [ i a t + £ t d^dr < f ||r A ||i 2(n;5 3) dr+11^^111,^3) 

Jo Jn Jo 

+ lk A (0) + C(^(0))||l W3) + /‘lie? IIL 2 (f2;5 3 ) d T 

Jo 

+ Hkt k)lll 2 (n ; s 3 ) + C(T ), (4.13) 

where the constant C{T ) does not depend on A and v > 0 is any positive constant. 
Theorem 14.11 and (14.111) imply that the first, second, third and fourth term on the right 
hand-side of (14.131) are bounded. Moreover, the assumption VM(dev(T(0))) G L 2 (fl; <S| ev ) 
and the convexity of the function M yield that 


/ M A (dev(T(0)))dx < / M(dev(T(0))) da: < / VM(dev(T(0))) dev(T(0)) dx .(4.14) 

Jo, J q J n 

Combining the inequalities (14. lip - (14. 141) the equality (14.9p becomes 


M x (dev(T x (t))) dx + 



o Jn 


C-'Tffidxd t+- / C(e t A (t))k(t)dx 


< C(T) + u\\e x {t) II l 2 (n ; s 3 ) 


(4.15) 


where the constant C(T) does not depend on A. Choosing v > 0 sufficiently small we 
complete the proof. □ 

Remark: The sequences {# a }a>o is bounded in the space L 2 (0, T; R 1 (0; M)) and the 

sequence {# a }a>o is bounded in the space L 2 (0, T; L 2 (h2; R)), hence it contains a subse¬ 
quence (again denoted using the superscript A) such that 9 X —$■ 9 a.e. in x (0, T). The 
continuity of / and Ti yield that 

e 

f(Ti(0 x + 9)) -f(Ti(6 + 9)) ^0 a.e. in ft x (0 ,T) 

and \ f (TA6 X + 9)) - f(Ti(9 +9)) \ is bounded independently of A. From the dominated 
Lebesgue theorem we conclude that for all q > 1 

f(Ti_(9 x + 9)) -f(Ti(9 + 0)) ^0 in L 9 (0, T; L q (Q.] M)). 
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Theorem 14.11 and I4~2l yield that the sequence of stresses {T a }a>o is bounded in H 1 ( 0, T; L 2 (f2; S 3 )). 
However, this information is not enough to pass to the limit in the system (13. ip as A tends 
to zero. We need to improve the convergence of the sequence {T A } A>0 . 

Theorem 4.3 (Strong convergence of stresses) 

Let us assume that the given data satisfy all requirements of Theorem [TH . Then, 


CT 1 (T A - T P ){T X - T p ) dx —> 0 


in 


in 


for X,/r —y 0 + uniformly on bounded time intervals. 

Proof: Calculate the time derivative 

j t (^ J C(e x - £ p - (e p,A - £ p ’ p )) (e A - - (e p ’ A - £ p ’ p )) dx) = 

' C(e A - £ p - (e p ' x - £ p ’ p )) (e x - ef - (e p f x - e p ’ p )) dx = 

C(e x -£ p - (e p ’ A - e p ’ p )) ( e x - e p ) dx 

G A (dev(T A )) - G p (dev(T p ))) (dev(T A ) - dev^")) dx 

= f (a x -a p + (C(e x -£?)(e x -£?))){e x -£?)dx- ! C(e A - ef)( £ t - e%) dx 
Jo Jn 

- [ (f(Ti_(9 x + 9)) - f(Ti(9 p + 9)))(dwu x -divuf)dx 
./o 

- J (^G x (dev(T x )) — G' p (dev(T At )) j ^dev(T A ) — dev(T^)j dx . (4.16) 

Using the fact that the given data for two approximation steps are equal and integrating 
with respect to time, we conclude that 


1 

2 


C(e A -£ p - (e p ’ A - £ p ’ p )) (e A - £ p - (e p ’ A - e p ’ M )) dx 



' 0 JO 


0 .to 



+ / C(e}-e?)(e}-e?)dxdT<D / |/(r,(0 A + 9)) -/(T, (0" + 0))| 2 didT 



0 «/o 


G A (dev(T A )) - G p (dev(T p )) (dev(T A ) - dev(T^)) dx dr , 


(4.17) 


where D does not depend on A, v. The remark before Theorem 14.31 and standard methods 
for maximal monotone operators finish the proof (cf. [2], [4] and urn- □ 


Proof of Theorem 11.2 


From the definition of the Yosida approximation we 
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obtain 


/ l|X(dev(T A )) - devfrjUj^n^dT < 

f IIX(dev(T A )) - dev(T A )f 12(!i;5lv) dr + f || dev(T A ) - dev(T)f 12(!i;5lv) dr 

J 0 J 0 

= X f ||VM A (dev(T A ))|| 2 2(a . 5 3 e jdr + f || dev(T A ) - dev(T)|| 2 2(n . 5 3 eJ dr . (4.18) 

J 0 J 0 

Theorem 14.21 implies that the sequence {VM A (dev(T A ))} A>0 is bounded in L 2 (0,T; L 2 (D; S^ ev )) 
and from Theorem 14.31 we conclude that the sequence {T A } A>0 is a Cauchy sequence in the 
space L°°(0, T; L 2 (D; <S 3 )). Passing to the limit in (j4.18j) with A —» 0 + we have that the 
sequence J A (dev(T A )) —> dev(T) in L 2 (0,T; L 2 (D; <S| ev )), hence it contains a subsequence 
such that J A (dev(T A (x, £))) —> dev(T(x,f)) for almost all (x,f) G (0, T) x Q. We know 
that 


£ p t ’\x,t) = VM(j A (dev(T A (x,f)))) = | J A (dev(T A (a;,t)))r 1 J A (dev(T A (x,t))), 
therefore 

ef’ A (x, t) —» | dev(T(x, f))| p_1 dev(T(x, t)) for a. a. (x, t) G (0, T) x D . (4.19) 

Now, multiplying equation (13.Ilk by v G i7 1 (0;M) and integrating with respect to D we 
obtain 

[ 6 x vdx + [ V6> A Vxdx + [ f(Ti(6 x + d))divu x vdx = I Ti(dev(T x ) e p ’ x ) v dx . 
Jn Jn Jn e Jn e 

The information f)4. 19)1 implies that 

Ti ( dev(T A (x, t)) £ p ’ X (x, t)) — Ti (| dev(T(x, t)) | p+1 ) —> 0 for a. a. (x, t) G (0, T) x 12 
and moreover 

|Ti(dev(T A (x,t))e£’ A (x,t)) - Ti (| dev(T(x, t))| p+1 ) | <^. 

From the Dominated Lebesgue theorem we conclude that 

Ti(iev(T x )E r i '\x,t)) — Ti(|dev(T)|» +1 ) ^0 in L 2 (0,T; L 2 (Si;S| ev )). 

Let us recall that (remark after Theorem 14.211 

f(Ti(6 x +~e)) -f(Ti(6 + 0)) ^0 in L 9 (0, T; L 9 (D; M)) 

for all q > 1, hence the function 

f(Ti(d x + 9))v -G f(Ti(9 + Q))v in L 2 (0, T; L 2 (D; R)) , 

where v G iL(f2;R). Finally, the energy estimates (Theorem 14.11 + Theorem 14.2[1 give us 
the following information: The sequence {w A ,£ p,A } is bounded in iL 1 (0, T; iL 1 (A2; M 3 ) x 
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L 2 (D; «S| ev )). The above informations are enough to pass to the limit in the Yosida ap¬ 
proximation and get the solution in the sense of Definition 11.11 □ 

Remark: The regularity of u t G L 2 (0, T; YY 1 (0; R 3 )) implies that the function 

6 t G L 2 (0, T; L 2 (D; R)). From the regularity theory for the parabolic equations we con¬ 
clude that the heat equation in (I1.6IU is satisfied for almost all (x, i)eOx (0, T). Hence 
the solution from Definition 11.11 is the L 2 - strong solution - for the definition we refer to 

0 - 
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